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SUMMABY 


!IIhe flow field in the close vicinity of -aie propeller in static 

V 

condition o has been analysed both experimentally and Iheoretioally. 
fhe static thrust of the propeller hf^ been measured by the wake survey 
method as well as by direct mechanical me'Uiods. 

A complete solution has been obtained for all the three components 
of the induced velocity close to the propeller for an arbitrary harmonic 

variation of vortioity as fxmotion of propeller radius. Eie Fourier constants 
have been determined from the circulation variation obtained from experimental 
data for the propeller* These constants have been used to estimate the 
coefficient of thrust There is a close agreement between the predicted 

values and the one obtaimd experimentally* 
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CHftPHER 1 


Introdnetlon 


35ie characteristics of paxipellers in static condition becosae 
important to oalenlate the nosing over moments and in estimation of take off 
amns* It is, therefore, essential to obtain accurate static tharost and torque 
of a fall scale propeller with zero forward velocity. In a wind tunnel it is 
quite difficult to achieve zero advance ratio condition due to the fact that 
the propeller itself oaieates a considerable air velocity through the tunnel 
and is also subjected to the constriction effects of the solid boundaries of 
the wind tunnel. 


Diehl 932) has collected a great deal of static Ihrust data 
mostly extrapolated from wind tunnel tests for a large group of propellers. 

His calculations of static thrust are based on constants derived from 
Durand and Lasley’s'^ '(1920) tests on wooden propellers. Diehl ^^'(1932) 
has used these 3:results to esthsate the take off distances of airplanes. 

Some useful static thrust data is also given "by Louis 1937 )• Unfortunately 
the material given in these references is not sufficient for a detailed 
estimation of propeller characteristics « Driggs^^\l938) has developed a 
method by means of idiich test data for a propeller obtained in a wind tunnel 
can easily be converted to aerofoil coefficients representing the propeller 
based on the assmption that the differential thrust and differential torq^ts 

dCSif-J 



/ dca \ 

respectively when plotted against CO 


produce curves idiich are sesaiellipses. This simplifies the expression to 


the extent that the total thrust and the total torque coefficients are 
merely the product of and the derivatives / Jand J 

at S5 ■ 0*5 • station representing ihe idiole propeller is thus found 
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to be CO 




» 0.5 or CO 

rs 


0*707 and is for oonTenienoet 


taken as the station at 0*7H. Conversely the thrust and torque 

/cIGt \ /dCa N 

gradients|-^^^y O. 7 R station may be easily obtained if the 

total ^urust and torque ooeffioients C^ and are known from propeller 
tests. However f this theoiry limits itself to propellers without bodies 
inspite of the fact titiat this analysis has been used fairly successfully 
for wind tunnel tests of most of the propellers provided that the hub and 
oylindrioal blade shanks are covered by a spinner* Montgomery and 
Hafner^^^(l957) have oasnried out a very satisfactory analysis of the 
static thrust of lifting airscrews both experimentally and theoretically* 
fhe experimental results obtained by them agree well with theoi^etioal 
calculations* fhe theoretical aspect has been analysed on the basis of 
vortex theory with the assumptions 1±tat the ma'&ematical model has infinite 
number of blades with small inflow angles neglecting rotational and radial 
oo^onents of induced velocity and slip stream contraction* She main 
feature of tiiis work is that the solidily as an independent parameter has 
been eliminated* At the same time this analysis is valid only for an 
evaluation of constant pitch or constant chord rotOTO. 


SoB» research a^ncies in the late fifties started focusing 
Interest on the static thrust of propeller especially after the subject 
of propeller driven VTOL alTOraft started to receive a considerable attention. 
0oward^^^(l957) gives a general method for establishing boundary conditions 
irixereby a efficiency under both static and forward flight conditions 
may be realized. Coward' '^(1959) also has given some data 
which yield a rapid} simple and sufficiently accurate estimation of static 
tiirast of propellers of conventional design for most preliminary analysis. 
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!Ehe purpose of present analysis is to study the flow field in tiie 
elose vifteinity of a propeller using the meastirei&ents of the pressure and 
velooi'ty distributions behind a Eartzell metallic propeller coupled to a 
lyooming 250 HP engine mounted on a four place aircraft Kanpur II. Hie 
Tortex theory is used to obtain the circulation distribution. Static ttirust 
measurements have been made by tethering the aircraft. 

Hie flow field behind the propeller as such is of very oomplex 
nature* Hiere is a eonslderable difference between the flow in static 
condition and the flow that exists when the propeller is moving or when 
the propeller is placed in a wind tunnel alrstream. In the static 
condition the propeller tips are woricing in a stzong vortex that tends 
to reduce the effective aspect ratio of the propeller. Based on tiie 
vortex theory concept a great deal of work has been done on propellers 
and various correction factors have been determined. Mostly this analysis 
has been confined to the Betz '^'(1919) optimum condition for lightly 
loaded propellers neglecting the slipstream contraction and the periodicity 
of flow. Hiis is equivalent to assuming '^lat the number of blades is 
large. PrandtX^^^(l919) bad developed a correction factor for a finite 
number of blades. Goldsteia^^^^(l929) has found out an expression for 
the circulation for lightly loaded propellers dealing with the conditions 
in the final wake. His basic assus^tion is that tiie velooi'^ field of 
the vortex s;^tem at a large distance behind the airscrew is equivalent 
to the potential field of the rigid heliooidal surfaces* Hiis solution 
when applied to the general problem is rigorously oorreot only for 
airscrews where the distribution of eiroulation along the blade is such 
that the flow in the wake is identical with the potential flow of siush 
a set of equidistant ooaxial heliooidal surfaces* It also assumes that 
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■ttx* path of the trailing vortices y asking up helicoidal sheets coincides 
with the resultant velocity. Theodorsen' ''(1944) has shown that the 
vortex theory can be extended to heavily loaded propellers by appljdng 
a correction factor to the slips tTOsa contraction^ 

In the present paper effort has been made to determine the 
velocity distribution in the close vicinity of the propeller by, assuming 
the variation of circulation along the radius of the propeller in the form 
of a Fourier series. The Fourier coefficients have been computed from 
experimental data. 

In Chapter 2 of the present thesis a mathematical formulation 
of ihe static thrust problem of a propeller is described. A mathematical 
model has been developed based on the usual assximptions made in the past 
with some modifications. In Chapter 3 '^is mathematical model is solved. 
Chapter 4 has been devoted to study the proHem experimentally for 
correlating the results with Ihe theory developed. An expression has 
been fotind relating circulation with the rise in presstire behind the 
propeller* !lhe experimental results obtained have been discussed in 
Chapter 5 while the conclusions of the present work have been given in 
Chapter 6 together with sii^estions for further investigations. 


CEAPaJEE 2 


inalytloal Foimulation of Problem 

She flow field behind a propeller is of a hi^ly complex nature. 

It is difficult to fozmialate a complete mathematical model and obtain its 
solution^ 4^^,.ts entirety taking into account this complicated flow field. 
Satisfactory resets have been obtained for any design problem by resorting 
to various assumptions. 

Since the propeller has no forward motion under static conditions 
and further since the axial velocity of the air is small compared to its 
rotational velociiy ihe vortex filament shed from the propeller blade move in 
helical paths whose pitch angles are nearly zero. We may, Idierefore, assume 
that the wake of a static propeller is composed of an infinite number of 
concentric circular vortex rings of constant strength forming continuous 
semiinfittitely long vortex cylinder (Pig. 2.1.1). 5he case of helicopter 
rotors under hovering conditions away from the ground effect is closely 
similar to the case of a static propeller. This assumption is also 
equivalent to assuming an infinite number of blades fosmiing the actuator 
disc. However, when the propeller is in forwas^d motion this assumption 
may not be strictly valid especially for the usual operation range of 
of propellers. Close behind the propellers there exists a strong axial 
variation of the helix angle of the trailing vortex filam^its. It is 
further assumed that the contraction in the slipstream is negligible. 

We shall assume that the circulation be constant in Ihe azimuthal direction. 
!Ehus ensuring axial symmetry as is generally valid in the case of airplane 
propellers. However, this assumption is not valid in the case of helicopter 
rotors with forward velocity and propellers in flight moving in pitching 
or yawing attitude. In these oases it is necessary to take into account 
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both jE»dial and azinruthal variation of circulation. Drees ^^^^(1949) has 
found an expression for the induced velocity of forward moving lifting rotors 
accounting for nonunifoxm azimuthal vortioity variation assuming however 
constant vortioilgr along the radius and infinite number of blades. Ihis 
assumption is not significantly different from the actual case of static 
propelleors becanse of ttie low wake velocities. 

As mentioned above the slipstream of the static propeller may 
be regaairded as an infinite series of coaxial circular vortex sheets 
(fig. 2*1 .2) parallel to the plane of the actuator disk. !^e vorticity 
of -tile sheet laay be approximated by a series of closely spaced concentric 
vortex rings. 

We consider first the case of a propeller having uniform 
circulation along tiie radius. This will produce a single helical vortex 
lying' along a oircvilar cylinder (Pig.2.1.5). In the far wake the induced 
velocity witiain the cylinder is uniform and there is no velocity induced 
at all outside the cylinder. In the case of varying circulation for an 
infinitesimal annulus of the actuator disc tiie irake is composed of two 
ooncentrio oiro^llar cylinders having eqtial and apposite vorticity and 
separated lydn*" • ®he dowmrash between the two cylinders is uniform and 
proportional to the annulus loading. Everywhere else the downwaah is zero. 

If the disk loading is made up of a continuous distribution of such annuli, 
ihe induced velooily dioreotly downstream of eaujh annulus is thus deteimned 
only that annulus loading and is not affected by the loading of any 
other annulus. Hence the radial variation of the induced velocity is 
proportional to radisil variation of the disk loading^ '(l956). 

Because of the variation of eireuLlation along the radius, 
trailing vortices spring from every point of trailing edge, thus forming 
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a Tortex sheet of helieoid shape, the zaimher of sheets heii]^ e^ual to the 

anabex of blades (fig.2*1*4)> The strength of each trailing vortex is 

determined from the ciroulation around the blade to satisfy Stoke 's relation 

of eontinuily of vorticity. Ifcj/^ is the increase of circulation between 

points )l~ and dt*" of the blade, the strength of the helical vortex 

springing froa this blade element will be— d/^. Therefore the total 

circulation/^ at any radius may be considered as the sum of all the trailing 

vortices idiich leave the blade outside that radius (Fig. 2.1 •3). As assxuaed 

by Heyson' ^'(1959) the close spacing of these vortices converts the 

slipstiream into a cylinder having a sheet of oontimous vorticity. Hence 

the calculated induced velocities will be only the mean or time averaged 

velocities. The axial component of the vorticity idiich is important in case 

of helicopter rotor blades with forward motion is assumed to be negligible. 

■b 

This is equivalant to sating that the rotor tip speed is infinite so that 

A 

the circulation is proportional to the local di^ loading^ ^■'(1960). Farther 


the blade bound vortex is neglected since its average affect in the ootirse 


of one revolution is zero. 


^e essential part of Ihe present analysis which is for a 
partioTilar propeller condition, is that the wake vorticity variation along 
the radius is measured and incorporated in the analysis. As already stated 
the wake cylinder will be filled with radial vorticity of a strength 
proportional to radial derivative of 1iie ciroulation. The equations which 
follow include all the three components of the induced velocity in the close 
vij^einity of the propeller. 



CHAPTER 3 


3*1 Method of Solution 

Let be the circulation per xinit length in the axial direction 
and it is expressed as 


UjCK) 


(5*1.1) 


stoere is the constant part of ciroTilation and X is a function of 
The elemental circulation per unit length along the coordi33ate % will 
be given by 



( 3 * 1 . 2 ) 


In normalised form with respect to and the total radius of 
propeller we have 



(3*1*5) 


Consider an element of a circular vortex ring of radius / at 
a vectorial distance . Let ■toe width of toe ring be J?. * 
circulation strength of toe ring can be written as — — • ■^s the 

dl 

oii^OBlatioa aloBg Hxe radius is also charigiug the strength of elemental 
circulation will be given ^ ^ ^ hi^er order terms 

of infinitesimals have been neglected* 


In toe first instance we consider only toe constant part of toe 

dTo . 


circulation i.e. the elemental ring with constant vortioi-ty 
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©le velocity potential <3> doe to such an element of vortex is 

. . (16) 

given hy 

^ = ±\[ ^ olA (5.1.4) 

^ atiJJ d)! 

'Qaantiiy £^otA measures 'Uie elementary solid angle, 

fherefore the velooiiy potential at point P can be written as 

where is the solid angle. 


<$> ZL i 

47r 




Adi 

dj! ' 


OT 


<5> _ ^_L[/8 CW -/3CJ?0] 


(5.1.5) 


We consider the point P at the plane of rotor. ©len for the entire 

cylinder (I, - O shid — c»0 


fherefore » 0 stnd jh 

Substituting these values in (5.1.5) we have 


cb — 1 i 

^ ^ djt 

80 that 


Wo r: ^ 

ct^ ^ 


or in nondimensionalised form we have eq» (5«1.6) as 

UJq — _ jL To 
COR ^ 


( 5 . 1 . 6 ) 


(5.1.7) 



Puttln# Xj-TtCO lire get 


ciS - dU) ~ 

K jr 

fox the limits when 

t ■=: O Zo ~ O IC - o 

tz R ^ '=: i_ *)(; " 7 ^ 


(5.1.8) 


(5.1.9) 


to 

We shall assume the oiroulatioVhe irepresented hy the Fourier series 

r' - /o^ C?c; r: [Ao 4- Av ■+ A^z. - --- 


with the Feujjier constants laormalised with respect to the circulation 

dT ^ TAo '- t- A\Sty\%~^ Ax^sr\2X-^ AvnSih mvl (5.1,10) 

d? di? ^ 


Bow 

6 (dr \ rr fL dr \ ^ 

j-r 

ct r r d (/y ) ^ 

“'idl)'' d(-A) J 

^ J7 d_ (^\ (5.1.11) 

P dG 5 ^ 

Substituting fordto ®<l* ( 5 * 1 * 8 ) in equation (5.I.II) we have 



XJsiBg eq.(3*1*2) we filially obtain 


d f cW s 

dt ^ 


V 


/? clC 

d^' 




1Wiere,.|'(^ %^ denotes derivative of Cy:>) with respect to ^ 

~ ^ ^ ^ “ 7T _S ,dZ? (X^Cosm^ 

^djl' / R2- Z_ ^ 


liters 


Ovyi — nn 


(5.1.15) 


Sow consider the effect of total circulation, fhe velocity 

induced by such a vortex riaig at a point P is given by 

( 17 ) 

Biot-Savart's law as' ' 



Rpom Fig. (2.1*1) we have 


<0 zr 1 (^V-cos'\j;j H- ^ L~A) 


therefore 


d? JC Co^jd'Vp 


(5.1.15) 


n- 


JL C y- Cos'V -><; -f jc Y-Sm-W'-y ) + -^) 


(5*1.16) 



12 


Su^atituting the above in eq. ( 5 * 1 * 14 ) we have for the total velocity 
induced by the irake 

i 1 


rrf 




47T 



_ Y'SrnW 

Y* cos 14) 

0 

y-cosiii -X 

Ysmw-Y 



(5.1,17) 


/ ^(-fcosyj ^ ^/2 


JJJ 

doe 


5*2 Galonlation of Component UJ of Induced Velocity 

Prom eq, (5*1 •I?) we obtain component of induced velocity as 


LU - 


L 

4rT 


an 1^ 

^ d /dr 


ofr" 




j j j 

000 


L^L. 

x/'cos U? +Y • dt\i; 

[ y*^ 4- „ 2 Cx sir» 


(5.2.1) 


Aftea? imtegratisg with respect to X suhstitutiiJg liidts and considering 

the flow field close to ttie propeller i.e, eq. (5.2.1 ) yields (Appendix l) 


UJ 


47T 


J 


^ " A (^\ 

olft V el^ ^ 


'o -^G 


os V -VY ^ dy# olty 

Y'^4-Y'Z. V - x7 T 


(3.2.2) 



fe aondiiaeasionalise the eq.(3.2.2) with respect to R and substitute the 
waives from eqs. (3«1*8 and 3*1*12) 



£ — -4- ^€o c 03 -t TpS5 SmV J (3*2*3) 

'25^ 4 f ^ 7p^ 4 - ;2 fto coS'M^* — Zp^ Sm^ 

On integrating with respect to '\|/ aztd substituting the limits (AppeMix II) 
eq» (3*2.3) yields 



Substituting for 


cL (^\ 


TT 

tU — .. ( y, 
4P. 



m^i 


fTOm eq. (3.1 *13) 


"m 


7^ 

Cos dX 

o 


or 

lu 

Ct>R 


i ^ 


m^o o 

i(v^ 


r 

Cos m% d ;)(; 


or 


(3.2.4) 


(3*2.5) 


(3.2*6) 



(3.2.7) 


This indicates that the incremental iwrmal component of induced velocity 
along the radius is zsto. 
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5»3 Oal0ttlatlon of X CoBtponent U of Indaoed Yeloei'fey 

From eq. (3«1*17) *e obtain the X coa^nent of induced velocity 


2 . 7 r 


U n — J- 


7 T ^ [^) 

df Velj^/ 

0 ^ o 


[ V"Co 314 / jj-tl) 


(5*5.1) 


f'^4- X ^ “Y ^ ~ ^ I*" cosip 4 - Y sTi r>ip ^-+‘2/2. 

As in the previous case ve substitute ahd integrate with respect 

to ^ and substitute the limits the above equation yields (Appendix III) 


fr^ r 

U— - ^ \ L V- Cos ip J elf- clip 

4 ‘n^ df \ df +72. C2>c/'cosiii ^-aV/s'in'H^ ^/2 


We nondimensionalise with inspect to R and obtain 


d / dr \ [‘So Go s ip] d^-clXf 


£ f 

22/3 


"^^^2 © Cos^' 

< ■« *=« 


(5*5*5) 


Proa appendijc IT integration with respect to'\|J yields 


\j - R ^ f d f ci_r N^'^dSo 
4J^ dC6 orfV 


( 5 . 5 . 4 ) 


Substituting: the value of ( — p- ^ y from (eq.5*1.l4) 

nhSS a C . 


d© ^ djZ' 
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u-=- — -E- s 

clil^ C¥%9'^) 



Ytsz^ nr 

^ (3-m. %^os m 


m-i 

-^0 

Fiob appesadix V ve have 

0 r: 

47 d;?' 

Yyv-^oo 

nn^-L 

or 

41TI Ijij;/ 


(5*3.5) 


^=0'. J_ ^ 

dj!' W^-fTp^ / 




(3.3.6) 


5*4 Oalottlation of y Component 19 of Induced Teloeii 


(5.4.1) 


From eq. (3*1.17) we obtain -tiie y ooa^nent of induced 
velocity as 


b:= 

47r 


2.7T P oO 


_d_ f dT 

df Vctji 


O O 0 


__ _ CUi')']dt-dJi. du> 

[y- ■+ '5'^+Y^+Z^-2.r^C>‘C0SXp+y5,n 

It Mjr I» Been that e4.(5,4.t) dlffere fron eq.Cj.J.,) i„ „3pect 
®f only trigononetrio ter. SmV and Ccs^ In the nnnerator. 
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Hence aolution of eq*(3*4«l) can l^e irritten aa 





dr' 

d 


Co 


or 

U • i 
— 

2?r 


dr' y 

djl' 


la. 

md 




3»5 Caloidation of Throat 

The tangential component oif the reaultant yelooity 
element is GL>f^ (Fig. 3.5.1) 


clr n POj n cIy^ 



From appendix 71 we hare 

\ 

■n YY'3'«^ 

-r r X 

m-1 

or Head imensionalis iiag eq*(3*5*2) by dividing both sides by 


We have 


Or = JI 

2. 7r 

[Arr+ 

Ava- 

m 

1 





^ / 

f* 



msoo 

^ hf7s 

XTT 

inhere s :i I 

Ao'TT 


2'\* Aiyv.*] 

ZL' IVL J 




m-i 


(5-4.2) 


( 5 . 4.5 


at the blade 


( 3 . 5 . 1 ) 


( 3 . 5 . 2 ) 


( 5 . 5 . 5 ) 
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She auiaber of circular vortex elements formed per unit time at 

The distance between the successive vortices is 

( 5 . 5 . 4 ) 

Since angle <:p is quite small as already defined we have (Fig,5.5*l) 


the rotor disk is 


:2.7V 


ci — 

a j- S.n <p 


<d r; tUo qTT tUo 

b b 


( 5 . 5 . 5 ) 


dS . 



KHy) 


alTUJo 


( 5 . 5 . 6 ) 


In nondimensional form we have 


dr’- _ T-'V/V'I kojK 

_ _ _ .fcx; _ U -I-CX) _ 


^TT UJo 


( 5 . 5 . 7 ) 


Prom eqs. (3.5.5) and (3*5.7) we have 


dn 

di' 


Gp J_ 

UJo ^ 




Substituti^ the value of from eq.(3.1.7) we have 

dX' 

( G-p \ 


Z s 


(5.5.8) 


(5.5.9) 


In rotor parameter the thrust is given by 


J? b r ClC d t 


J 

o 


(5.5.10) 



ia 


Assuming Ci_-z.(X<S^ we can write eq.(3.5.9) in the form 

T •=: ^ Oo bco^JoC e 

Since 

0 - q:)) Pig. (3.5.1) 

,2 


X 


X. QobCo^ 


j 


Co ~ CpJ> c dY 


In nondimensionalised form we have 

I 


Cr 


dob 

:i7r 


J 


C^Q ^ qp') d 03 d te? 


(5.5.11) 


Gch 

2.7^ 


( Q ^ ^ J-)c.' 'So^dSS 

COR 5 ; •' 




o 


(3.5.12) 


Prom eqs. (3.5.9) and (3. 5. 12) we have 

y \JJo \ 2- ^ 


or 


4Tr 


A 


tUe 4 


O _ -m_ _ 


LUe Ciob S 


) 


4Tr 






r 


O 


(z‘S3d^ 


or 


^ flobs' 




J 


C ^Co cJco 


gobs 

4“^ 


/d 


0c'53^dto 


-2, 


Where 


OUq 
Cor. 

4- *4 A ^ ^ 


(5.5.13) 


(5.5.14) 


~ o 
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Solving the quadratic we have 

X - ~»4 

Z 


The quantity on the ri^t hand side of tiie above equation has to be 
positive 

X can be determiiied by evaluatingy^ and ^ irtiere 


.1 

A - <hk£ic'dicit3 
- CfvJ 

O 


o __ do b s \ 0 c^’co^cL to 

^ — zlfTT 


( 5 . 5 . 16 ) 


(5.5.17) 



CHAPfEB 4 


4»1 Experimental Setup and Meaeureaenta 

A propeller produces a pressure rise across the propeller disk, 
fhis may he related to the circulation distribution over the disk. We shall 
now obtain an expression for the pressure rise in terms of the circulation.^' 
Strip theory can be used to obtain an expression for the induced velooi'ty 
variation over a nonuniformly loaded propeller of arbitrary planform 
twist. Consider an annulus of radius and width off (Pig. 5.5.1) of the 
propeller disk. We shall assume the inflow single cp to be small so that 

5m cp n: -kl Ti cp “ ^ 


Mass of air flowiiag throu£^ Ihe annulus 

Ml:. 2.t\ Y~ 


Elemental thrust gIt- elemental lift 

Xhe total elemental lifter i>number of propeller blades 
— b-dL ~ 2^1^ CV|H-UJo) 2 Wo df 
or 


b — - 
d® 


XtX p Cyt -f- Wo) 2 U)o Sj 


(4*1.2) 


(4.1.3) 


In Pig. (4.1.1) station (l) is on -the far upstream of the propeller. 
Stations (2) and (5) are the positions shown close to the propeller disk 
upstream and downstream respectively. Station (4) is on ihe far dowzistream 
in the propeller wake, therefore we write 
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It, - f^o -+ ± />V,2- (4.1.^) 

^ — |^° ■+ ^ = (^O + -t ' 2 .W 0 J ^ (4-1-5) 


Dhe rise in total inressore through the disk is 


^P- ^ — Rtjl r’ 1 fc^tv^, -f 4 u;o^^ z puJc ( 4*1 *6) 


Henoe from eq.s. (4*1*?) a^d (4*1*6) we have 

^ 4 ^ ^ ARr 2 ir 


From vortex theory we have for tiie lift coefficient 


Ol n 

■JL C \/y ‘Co 
X 

So that 



dL 

g (‘55 


« — . H f\\ ! 

- -— f^r CO 


( 4 * 1 * 7 ) 


( 4 * 1 . 8 ) 


( 4 * 1 * 9 ) 


eq8*(4*1*7)t (4. 1.8), and (4.I.9) we have for the circulation at 
radius 


ZlRf 

P 00 b 


(4.1.10) 


Assuaing circulation and as function of 'cS have from eq.( 4 * 1 * 10 ) 



Co 

C^oJ ^ fb 
xn 

(4*1*11) 

We shall 

define a pressure coefficient Cp for liiis 

oase by dividing 

both sides of (4.1.11) by — pco so that 


Cp = 

4PtCS; _ rtsTj b 

COR^ 7X- 

(4*1*12) 

or 



— 

b PcS) _ b r^SJ 

~ 7r 

(4*1.15) 

where 

P - CJR^ 


since 

APt - &<jk 

( 4 * 1 . 14 ) 

- t Fmj 

^ fco^ v: r 

2, 

( 4 * 1 . 15 ) 

or 



Tcsj 


(4*1.16) 


oo^ 

wheire 

^ _ xrrfu<^ 

( 4 . 1 * 17 ) 


Xhis establishes a disrect relation between the oixculation and 
the pressure rise across the propeller dii^# 
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From eq* (4*1*7) we obtain for the elemental thrust d 7 over the 
element d f" 

cIt n A9j ZTT S3 d Sj 

or 

ch - 2.TrR^S3 /g, ghdto 

(4*1.18) 

Oliis can be Integrated over the radius to get the total thrust. 

4*2*1 Measuroment of Static lEhrust by Wake Survey 

A 7* diameter sietallio two bladed Hartzell propeller installed 
on a four place single engined aircraft Kanpur II has been used for this 
purpose. The propeller was conneeted to a Lycoming 250 HP engine (Pig. 4*2*1). 
Bffort was made to obtain blade profile specifications and propeller 
characteristics from the manufacturers well in advance but the same could 
not be piK>oured in time, fhe i^asured variation of the chord and blade 
angle has been plotted in (Pig. 4*2*2). Sixteen brass tubes of outer 
dia 0*069 iii* and 6 in. long were mounted uniformly spaced 3 in* apart 
on a wooden strut spaiming the oosplete radius of the propeller with four 
of the tubes projecting beyond the propeller disk* The rake was mounted 
on the top of the engine nacelle (Fig. 4*2*3) parallel to plane of the disk 
and at a distance of 3 * 0.5 in* behind the plane of rotation on the axis 
of "T! * Itoe to the safety requirements and installation difficulties a 
closer approach to the plane of rotation would eliminate the effects of 
tip vortex roll up in the outer wake for^^ — or the static propeller. 

n O “ 

After imsunting the Pitot rake on the aircraft the position of the tubes 
along the radius was measured. Fifteen tubes (2-1 6) on the rake were 
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oozmected to a maltimanoioeter vith fifteen glass tubes* fhe fii»t pitot 
iras oozmected. to isiciromazioseter* Azzother micaxzioanoaetex was cozmeoted to 
a tube fixed ozi the wizig quite far fi»m the propeller to observe the ohazsges 
in the djrziamlc head of the free stream* The danipizig inhezrezit Izz the 
eombizzatiozz of pressure imbesi aaizometer azzd the cozmectizzg tubiEig was 
sufficiently large to allow the assuz^tion that the manometer zreadizzgs 
represent the steady time averaged pz?esszu^s* 

The multimanometer and the micromanoraeters wezze kept in the coctpit 
(Fig. 4*2.4)* After levelllzig it latezzally the multimanometer was szispeztded 
with the help of a silk thread to ezzsure automatic verticall'^ while *^6 
pTOpeller was ronnizzg* The two miczomaziometers were fixed in the cockpit 
and levelled before startizig the engizze* Distilled water was used in the 
manometers to indicate the pressure chazzges* 

The tests were conducted by zruzinizig the propeller at ei^t diffeirent 
rotational speeds (BPM 750, 1000, 1250, I5OO, 1750, 2000, 2250 and 24OO) with 
the blade pitch locked in the fine position. 

4«2«2 Meastizyement of Chord and Blade Angle 

An improvised method has been used to find the blade azigle and 
the blade ohozrd at the various bl€ide statioz^* The aircraft was levelled 
up in lateral axis only* A big metallic plate (size 8*x4'xi”) was kept 
under the propeller on the ground and properly levelled* With the help 
of a pluffibob the propeller nose centre was fixed* This was fizz^ly used 
to place the propeller in a horizontal position* Projection of the farthest 
point on the plate from both the propeller blade tips adjudged by observation 
fixed the referezzce lizie on the plate placed on the ground. The propeller 
was then positioned in the lozigitudizzal plane by setting the lower blade tip 



at th« ndninuia distance from the plate on the ground. At equ^ intervals 
projections from ttie plate were dropped on the leadiiag and the trailing 
edges which in turn decided the farthest points on the aerofoil section of 
the blade* Finally the distance between these two points was measured with 
the help of vernier calipers giving the chord length at the section. 1316 lines 
joining the two points on the ground plane corresponding to the blade chords 
weire drawn and the angles with reference to the line joining the blade tip 
pTOjeotions in the horizontal position were measured. The variation of blade 
chord and blade angle has been shown in Fig. (4.2.2). 

4*5 Measurement of Static Thrust by Application of 
Miieot Mechanical Methods 

Two methods have been employed to measiire the static thrust direct* 

In the first case three dynamometer '^e spring balances of 200 kg. ^ 100 kg. 
and 100 kg. range were connected in parallel rigidly. The complete fitting 
with the help of ^ in. SWR was connected to Ihe tail wheel mounting (Fig.4*5«l) 
aikl a peg fixed in groxind ^not shown in the Pig.). propeller was run at 

Idle same rotational speeds as mentioned in Sec. (4*2.1) and the thrust was 
recorded as indicated by the dynamometer spring balances* To ensure the 
reliabili-^ the spring balances were calibireited before haial. 

In the second case the static thrust was measured with an SR4 strain 
gauge fixed on a calibrated (Pig. 4*5 *2) mild steel specimen. She calibration 
curve of Ihe strain gauge tension is shown in Pig. (4»5«3)» 

The thrust indicated by both the methods has been shown in 


Hs- (5.2*5) 
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4»4 Wake Smrrey by Pitch and Yawaeter 

A 2 in* diaiaeter five tube sphezre tjrpe combined pitch and yawmeter 
was used to detezmine the direction of the airstzream behind the propeller 
(Fig.4*4.l). The oombi^d pitch and yawmeter was made out of a 2 in. diameter 
spherical brass ball. On the circular cross section four holes were drilled 
for pitch and yaw measurement and the fifth hole was in the centre to observe 
the total head. This pitch and yawmeter was mounted on a strut having 
stireaffilined shape fixed behind the propeller. A glass tube was fixed on 
the rake close to the yawmeter to measure the local static head. 

3he measurements have been taken by running the propeller at 
five different speeds (BBC 750, 1000, 1250, 1500, 1750). The pitch and 
yawmeter was shifted at four different positions on the rake along the 
radius. The pressure diffemeces were observed on Ihe mioromanomters 
Installed in the airplane cabin. The ccHnbiiaed pitch and yawmeter was calibrated 
in the 5**2’ low speed wind tunnel (Pig.4.4*2) of the department. The 
oalihration curves of the yawmeter are shown in Pig. (4.4»3)* 



CHAPm 5 


Msottssion of Results 


5*1 Flow Field Behind the ProTseller in its Close Yioinlty 

It will be obserred fzpm eq« (3<>2*7) that the increioental normal 
component CCU9 of 1'ho induced velociigr along the radiue of the propeller disk 
is zero* Zhis indicates that the normal component of the induced 
velocity caused by the constant part of the circulation remains constant 

\ c|P 

in the wake along the propeller so that LU© zi. 3*1.6). 

Therefore it may be inferred that the normal oonqjonent of induced 

velocity in the propeller wake does not depend upon the variation in 
circulation along the blade radius. Ihis condition will not occur in the 
case of propellers in pitch or yaw and helicopter rotors in forward flight. 
If in eq. (3.1*5) choose a point P in the far wake the total solid 
angle subtended at p will be^TT so that the normal component of the 

ciP 

velocily induced will be "Zi 2. tUo • ^ quite close 

agreement with toe momentum theory of toe propeller. As seen from 
Fig. (5.1.1) the predicted static thrust is in close agreement with the 
measured static thrust for toe wake assumed in Chapter 2. 

Eqs. (5*3.6) and (3*4*5) show that the vertical and lateral 
ooiiq>oxients lA and\9 respectively are symmetrical about the y(, andV 
axis as has been assumed. Moreover both these ooa^onents of the induced 
veloci'ty do not play any part in the estimation of static thsrust* 

From an estimation of snd IjUo fl* 4® possible to calculate 

the total velocily induced in the wake and hence the elemental pressure 




dist3?il3atloa aloai^ the radius of the jojopeller la statlo oondltlon eaa be 
detemlned* 


5*2 Estliaatlon of fiEopeller Statlo fhrust 

fhe paresstare rise behlM the propeller la Its close vl^olalty 
obtaiaed from the wake survey has beea shoua la Pig. (5.2<»1). She 
measured data are subject to the follovlog oorreotloast 

(a) Maaometer Reservoir Correotloa 

(Che inoltlBianometer used has flfteea 0*145 la* diameter bore glass 
tubes ooaaeoted to the two 0*5 la* diameter bore storage alumlalum tubes. 
I«t l^iy )^3L — Kis-be the change la the level of the glass tube aad H be 
the change la the liquid level of the storage tube. Shea we have 



(5.2.1) 


( 5 * 2 . 2 ) 


fherefo3?e the actual rise la pressure la one tube will be 


h 4 - H 


(b) Oorreotlon for Slipstream Gontraotlon and Rotation 

fhe rake should be lastsilled on tiie aircraft behind the propeller 

be 

in such a way that the probes should /as close as possible to the blade trailing 

edge* this will avoid the effects of slipstream contraction aad rotation 

especially for V — n as under such condition the propeller blade tips 
riD ^ ^ 

produce a very strong vortex field* In -Uie present experimental analysis 
due to the installation difficulties the distance between the blade trailing 
edge and the probes could not be reduced below 3” i 0.5”* ffiiis is quite 
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large to cause oo^idei^able inaccuracy in the static thrust# HoweTer 

(op ^ 

Look' '(1939) has accounted for the slipstreaia oontraotion. A method for 
the estimation of the correction factor for slipstream rotation has been 
given by Look' ^(1952)» Because of Uie nonavailability of these two reports 
it has not been possible to correct for rotation and contraction of the 
slipstream. 


(0) Flow Inclination to the Pitot Comb 

The total head readings of the tubes of the pitot comb have to be 
corrected for the local measxired flow inclinations. ®iis requires a 
calibration of the total head tube with yaw. For want of time this has 
not been done. The correction will be incorporated soon after Uie 
calibration has been made. 


Pig. (5.2.2) shows the circulation variation over the blade 
(eq. 4.1.16) corrected for manometer reservoir level variations only. 

rr 

An arithmetical mean dimensionless circulation — has also been shown 
in Pig. (5.2.2). The shaded portion in this Pig. (5.2.2) indicates the 
blade tip losses due to the rolling up of the blade tip vortices and the 
wake contraction. In normal cases when the propeller has some advance ratio 

'::±. O the circulation curve should intersect the ^ axis between 
T»D ^ 

■ 0*95 and 1. She present experimental result substantiates the 
foregoing statement that the tip vortices are stronger in case of propellers 
with static condition as compared to propellers in forward flight. 


For computational purpose the mean curve has been approximated 
at the root and the blade tip as shown. Due to the presence of the engine 
nacelle the first probe fixed at GO - 0.21 had too many nonperiodic fluctuations 
at all propeller speeds possibly due to the engine nacelle interference. 
Consequently the correct behaviotir of air stream at this station could not 
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“be ottained, Again piofces situated beyond ihe propeller disk radius indicate 
a finite pressure rise possibly caused by the blade tip flow. However, Idhis 
appairent residual circulation has also been neglected for computational 
purposes. She approximated mean curve has been xised to compute the Fourier 
constants Abusing the least mean square fit. !Hie computed Foxirier 
constants (Appendix VII ) have been used to estimate theoretically tiie 
coefficient of 1iii?ust(GT) and normal component of the induced 

velocity using eqs. (5.5.9)f (5»5.15)» (5*5.16) and (5*5.17)^ and ^ 
have been obtained by numerical integration using Simpson's rule taking 
the measured values of(c') and(^^J) from Fig. (4.2*5). 

The static thrust estimated and those measured at ei^t different 
rotational speeds of the propeller using the three different methods of 
chapter (4.2*1 ) and (4*3) has been plotted against the rotational speeds of 
propeller (Fig. 5*2.3) • Figi(5«2o4) indicates -fee variation of thrust with 
respect to square of angular velocitjij^l^J of the propeller. Approximate 
correction for the 20 kg. dead load of the spring balance complex has been 
applied in respect of static thmist obtained by spring balance. Ground 
friction of 15 kg. has been accounted for in case of static thrust obtained 
by the spring balance method and the strain gauge method. Satisfactory 
linearity has been obtained with all the three methods. 

Coefficient of thiu.st(tr3 calculated and theoretically predicted 
has been shown in Fig. (5.I .l ). As may be observed both the values are in 
close agreement. It may, liiezrefore, be stated that the present analysis 
gives a reliable estimate of the static thrust of propellers having varying 
chord and pitch distribution along the radius. 

5*5 Flow Direction Behind the Propeller 


It may be obsesjved from Pig. (5*3*1 ) that the air stream has a 
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tendency to aove inward indicating slipstream contraction. !Ehe contraction 
reduces as the BIW of the propeller increases and decreases as we go along 
the radius of the propeller from root to tip. !Ehis is quite consistant with 
the assumed wake in case of the propeller in static condition. Shis is quite 
distinct in Fig. idii<h indicates the inclination of slipstream both 

in yaw and pitch with respect to the plane of the propeller at warious speeds. 

It may be added that effort was made to measure bolh pitch and yaw 
exactly at the blade tip. fhe same has been found to be almost zero due to 
the slipstream contraction. However from Fig. (5*3*1 ) it is possible to 
estimate the rotational velocity in the slipstream of the propeller as shown 
in Pig. ( 5 . 5 * 2 ). 



GHAPEER 6 


6e1 Gonolusion 

A study of syrametary relations concerning the induced velocities in 
the plane of a nononifornily loaded propeller in static condition indicates 
i^tt 

1) * ®ie eoB^nent (iUio) of the Induced velocily in the axial direction 

of the propeller in the wake is constant along ihe radius and it 
does not depend upon the variation in circulation as a function of 
the propeller jradius. 

2) . The induced velocity components in X and V direction are 

symmetjcioal idiioh is oonsistant with the assumed wake. 

3 ) . The static thrust of Ihe propeller can qtiite accurately be computed 

in case the circulation variation as a function of the radius of 
the propeller be known. Conversely the pitch and chord of the 
propeller can be varied to obtain the desired static thrust. 

6*2 Suggestions 

It is possible to extend the provisions of the present study to 
the case of helicopter rotors in hovering condition as the flow field in 
the case of propellers in static condition and helicopter rotor blades in 
hovering have approximately similar characteristics. The experimental 
procedure with some modification can be utilised to study the characteristics 
of the wake and thrust developed by the rotoi^ while operating in-ground effect 
and out-of-ground effect. 



APPENDIX I 


Solution of Integral 5»2.1 

■j - ir^^ xY-cos^r di 

^ C^<^os\^ ^ y 5'm Yj-f J Vi 

We make tiie folloviz]^ substitutions 


( 1 ) 


A ^ 

a ^ 
b - 


2 

f -f xy^cos i|f -t- Y|i sm 


Y ^ ” Zy c x Co S If f V s l“n\jr J 


o 


( 2 ) 


c = i 


4 o c ~ 


Hence eq.(l) reduces to the form 


I - 




o 


(a 4 bi + Ci^JZl 2 . 


6i 


(3) 


Prom item 16? ireference ( 17 ) 


I = A[ 


ZCci -tbj 

q JoTbFTr 


oO 

] 

O 


(4) 


Substituting the limits t eq* 4 yields 


I 


A 


C%4^c + b)Ja 


(5) 
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APPENDIX I (Continued) 


Prom eq. (2) we have 



a 

— COS j-Y i^ sm\y ) 

^ “i" ^ XV'^V. Co S H" V S”l T) V|^) 


( 6 ) 



APPEiroix n 


Eyalaation of Integral ^.2.^ 

JO <r'^ 

I L — CO -t~ 5 to C OS <>y •+ p a> si n-qj- J elixir 
& CO ^2,%^ cos c!5 S' i r\ 


We make tiie following substitutions 

A = 

B = 

c = '>p €o 

0- •= f^-i- rp- 

b = ~ 2, i to 


Hence eq.(l) reduces to the form 
2!Tt 

J =: I A i- l3 CosM> -I- C S\Y^W 


Q “f bcosiij 4-/^ 


cl Y 


i^m item 2.558-2 reference 18 we obtain 




O 4- bcosiH 4^ 4 


BUce 

b"4>C-^ 


•4 - 


4 Qc ^ \ 
b^4/:^ 




Q 4 b cos 14? -f ^ Srr^ijj 

•4 



IPEEiroiX II (Oontiimed) 


O -4- bOoSY '4-/- Sf r>M3 

Vf ^ z 

7b^+X:"-a^^(Q-bJ) i^tt +^y^i+c^_o^ 

• OTr 4n ^ -i- /C^ , X ^ -7 


- i 


— loj (^a-t-ye-l^ ^ ) tuhen Q zb. 


C5) 


+ C Oi-bJ-h^ yf 


loh^-y^ O^nb^C^ 


In all the oases stated above integral ( 5 ) will be zero after substituting 
limits from 0 to^'i^ 3!herefore eq.(4) reduces to the form 


i-[ 




On substituting the limits we obtain 


] = Bb -t C/C -..p 

-h /z^ 


substitute values fTOm eq.(5) 


2 Cs'‘'<^\ V^cS^) 


-TT 



iH>EKDIX III 


Ihralnation of Integral 


1 = 


C 'f Cos yf ) 


II V'^ H- X^-V 2. X eos'ip -V Y S4r» V 


7e malce tbe folio wixig substitutions 

A = COS 'll/ 

z 

0 " A “f-x^ -ty^ — Cx cos 14 ? -+ y 

b- = o 
C =: i 

q = 4 ac 


Hence equ(l) reduces to the form 

Ag- clil 

'o Co + bJ! 

From item 175 reference 17 we have 


ly 

**f*k 


iq- 


2 cM -i za) -Q 

'VoThTfcgi o 


I = -4c 


to 


xa 






Substituting the values from eq.(2) we have 
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da 


^ Cos V 


yTy- x^-v‘y^_ 



^IV ' 5 >|D 


APPEHDIX IT 


Evalttatlon of Integral 5o3«3 


j* (^COS 'V ) 

^ “+ t 2.%W COS"^ -2T?€u S*<r»'Vp 


We siake the following transfoianations 



C COS Ip -t ip Sl'n -ip? ^ 'X 


iEherefore 
^13 

1 




/^y/ •+/^K 
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or from reference 19 we have 


i T\3.^ ^ I 

I = i- [cos-q^r X pT^C/^-^'" ] 

r):io 




where 
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PiC /0 - i(^ 3 /^"-i 3 
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APPENDIX IT (Continued) 


^ 

^ t Ai) - — (jLZ\p^ — \OS~J^^-i>') 


W0 substitute (4) iu eq^#(5) and integrate term by teimi 
2.11 


1 =j COSTjf^l +/aK + £ y. 
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“ zn 

I 2 , - cos''V ( 'gcos'M) ■+ Tp sm'M'J olTl-' 
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APPEKDIX 17 (Continaed) 


Like-wise terms of all hi^er order are zero. Therefore, 


I 




( 9 ) 


or 




1 / 




(10) 
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APOTDIX V 
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IPPEHDIX VI 
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APiENDIX VII 


Following axe the Fourier oonetants as computed. 



0.00 

0.214 

0.286 

0.315 

0.557 

0.428 

0.5 

0,57 

0.645 

0,71 

0,786 

0.86 

0.95 


1 .755956 
-0.178640 

0.122000 

0,385681 

0,200718 

0.479401 

-0.065506 

- 0,006008 

-0.448655 

-0.306022 

-0,462631 

- 0.255855 

-0.172175 

-0.137559 


1,00 
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